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Supcrsymmetry (SUSY) in quantum mechanics is extended from square-integrable states to those 
satisfying the outgoing-wave boundary condition, in a Klein-Gordon formulation. This boundary 
condition allows both the usual normal modes and quasinormal modes with complex eigenvalues u). 
The simple generalization leads to three features: the counting of eigenstates under SUSY becomes 
more systematic; the linear-space structure of outgoing waves (nontrivially different from the usual 
Hilbert space of square-integrable states) is preserved by SUSY; and multiple states at the same 
frequency (not allowed for normal modes) are also preserved. The existence or otherwise of SUSY 
partners is furthermore relevant to the question of inversion: are open systems uniquely determined 
by their complex outgoing-wave spectra? 
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I. INTRODUCTION 



A. Outline 



Supersymmetry (SUSY) in quantum field theory 
relates bosons to fermions. Its analog in quantum me- 
chanics is interesting in its own right, and relates two 
(typically one-dimensional) Hamiltonians H and H with 
the same spectrum of normal modes (NMs) , or the same 
spectrum apart from one state. The classical limit || 
relates a one-parameter family of Hamiltonians -ff (£)• 

This paper generalizes SUSY in quantum mechanics 
to open systems and in particular to their quasinormal 
modes (QNMs), which are eigenfunctions satisfying the 
outgoing-wave condition (OWC) at infinity. The gener- 
alization itself is straightforward, and we focus on the 
novel features that ensue. After a brief introduction to 
open systems and QNMs in the rest of this Section, this 
paper presents three main results. 

First, in the conventional discussion for the 
Schrodinger equation, spectrum preservation can be con- 
veniently expressed in terms of the number of NMs n(E) 
for H and n(E) for H at the same energy E . Namely 
that n{E) = n(E) for all real E except a privileged value 
Eq (usually chosen to be Eq — 0). At this value, the 
difference is given by the Witten index p] 



A(E ) = n(E ) - h(E ) 



(1.1) 



which can be +1, 0, or —1. For open systems with the 
OWC, time-reversal invariance is broken, and it is ap- 
propriate to consider a Klein-Gordon equation (KGE) 



time-independent equation, and distinguishing +lo from 
— lo (since reversing the two converts an outgoing to an 
incoming wave) || . Section [n] shows that the equality of 
spectrum for the KGE extends to the complex w-planc 
except at the two privileged frequencies ±0 = ±\fE^, 
namely 



n(uj) — n(bj) 



At ±0, one needs to consider 



A(±tt) =n(±n)-n(±n) 



(1.2) 



(1.3) 



instead — in effect replacing <9 t 



df or E 



the 



These are again +1, 0, or —1, but with A(S1) = — A(— SI) 
determined by the asymptotic behavior of the SUSY 
generator, to be defined below. In other words, under 
H i— > H, if a state is removed (added) at SI, then a state 
is added (removed) at —SI. This relationship, which also 
applies to conservative systems and NMs provided we 
take the KGE point of view, sharpens the information 
provided by (1.1). 

Secondly, SUSY preserves norms and inner products. 
However, for outgoing waves the usual norms and inner 
products are not useful. For example, an outgoing wave 
of frequency lj goes as exp[iw(|a;| — t)] at spatial infinity. 
With Im uj < for QNMs (see Section IB), the exponen- 
tial growth in |a;| renders the wavefunction not normal- 
izable in the usual sense. A generalized norm for QNMs 
was first introduced by Zeldovich |?J many years ago, and 
shown to be useful for time-independent perturbation 
theory (of the complex eigenvalues) An associated 
generalized inner product can also be defined |^|. The 
time-evolution operator turns out to be symmetric un- 
der this product (the analog of self-adjoint). Section III 
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shows that SUSY preserves these generalized norms and 
inner products — a pleasant surprise, since their con- 
struction is totally unrelated to SUSY. 

Thirdly, for non-conservative systems, there is no guar- 
antee that the Hamiltonian can be completely diagonal- 
ized; in general the best that one can do is to decompose 
it into Jordan blocks (JBs) [Q. Each block j, say of 
size Mj x Mj, is associated with an eigenfrequency cj 
with Mj = 1 being the usual case of a QNM. Section 
shows that SUSY preserves the JB structure: except for 
u>j = ±f2, a block of size Mj maps to a block also of size 



4(x) 



^+iuj\x\ 



OO 



(1.6) 



Beca use V is trivial outside J, the asymptotic conditions 
(1.6) can be stated at x = ±a instead: 



<j>{x) 



±a 



(1.7) 



The imposition of two boundary conditions in (1.7) 
forces the eigenvalues to be discrete, and these fall into 
two classes. First, there could be bound states or 
NMs |l4|; these must (from the Schrodinger point of 
Mj at the same frequency "ujj. In fact, if we generalize view) have E = uj 2 < 0, and hence u> is purely imaginary, 
the definition of n(u!j) to be the order Mj of the bl ock , 
then the relationship between n(u) and n{u>) [cf. ( |l.2| ) 
and (L3)] remains valid even for n,n > 1. 

Examples are given in Section [y] and a discussion is 
presented in Section VI , including a sketch of some issues 
for potentials with tails — which can lead to situations 
with negative n and n that are nevertheless accommo- 
dated in the same formalism. 



B. Quasinormal modes 

In open systems, waves are not confined, but can es- 
cape: acoustic waves from a musical instrument, elec- 
tromagnetic waves from a laser, and linearized gravita- 
tional waves from a Schwarzschild background (to infinity 
and into the horizon) . These systems are often described 
(e.g., in the case of gravitational waves H]) by the KGE 



[d 2 - &l + v(x)] #m) = o 



(1.4) 



or (e.g., in the case of optics |TJ]) by the closely related 
wave equation 



[ P {x)d 2 t -dl] 4>(x,t)=o 



(1.5) 



This paper deals only with the KGE, both because it 
is readily related to the Schrodinger equation in terms 
of which SUSY is usually formulated |^| , and also be- 
cause it (unlike the wave equation) admits NMs which 
are interesting in the pres ent discussion. 

Except for Section VI D , we shall assum e th at V(x) [or 
p(x) — 1 in the case of the wave equation (1.5)] has finite 
support on the interval I = [—a, a], which is natural for 
describing a system of limited extent, surrounded by a 
trivial medium such as vacuum. 

We assume that the loss is only due to the boundary 
conditions. In particular, the potential V is real. Ab- 
sorption may be described by a complex V, but causality 
then requires dispersion; the necessary generalization E3] 
will not be discussed here. 



Since bound-state wavefunctions vanish at infinity, (1.6) 
dictates that Im uj > 0. Second, there could be QNMs 
with complex eigenvalues uj 2 . Because these waves es- 
cape, 4> decreases, so Im ui < pBfl . Provided Re w ^ 0, 
they occur in pairs: u)_j = —Wj, as is readily shown by 
conjugating the defining equation and boundary condi- 
tions. Those with Re u> = need not be paired; these 
zero modes [ [l6| will be of particular importance below. 
Consider the potential shown by the broken line in Fig- 
ure 1(a); its NMs and QNMs are shown in Figure 1(b). In 
this example, there is one NM (triangle) and a sequence 
of QNMs (crosses), including a zero mode. In contrast 
to the Schrodinger formulation, the use of the KGE and 
the introduction of T allows NMs and QNMs to be dis- 
cussed together — and at least in this example the latter 
manifestly carry much richer information. 

Even though QNM eigenfunctions are not square- 
integrable and do not form a conventional Hilbert space, 
they are useful for analyzing outgoing waves. Impor- 
tantly, the complex QNM frequencies are often directly 
observable: e.g., the central frequency and width of an 
optical line observed from a laser cavity, or the rates of 
repetition and decay of a gravitational-wave signal that 
may within the next decade be detected by instruments 
such as LIGO |l7]]. In fact, the spectrum is often so rich 
that, under some broad conditions, namely that V(x) 
vanishes outside [—a, a] and has singularities at x = ±a, 
the eigenstates cf>j of T are complete for x £ I ]1S[ |, so 
that for t > the observed wave signal can be repre- 
sented as M 



4>(x,t) 



dj<l>j{x)e 



(1.8) 



Among the solutions of (1.4), we consider the space T 
of states satisfying the OWC. We leave the time-domain 
definition to Section III ; in the frequency domain, a func- 
tion is in r if: 



In cases where the ifij's are not complete, it may still 
be possible to characterize the remainder, which could 
be, for example, a power law in t for long times p0| , pl| . 
Moreover, when the eigenstates are complete, one can 
set up a formalism that closely parallels the conservative 
case (see Refs. P,p2[ and Section IIIB). One can even 
second-quantize using these eigenstates as a basis (e.g., 
to discuss thermal effects and atom-field interactions in 
an optical cavity) |^3| . These developments have been re- 
viewed |^] . We shall see that much of the mathematical 
structure is preserved under SUSY. 



2 



II. FORMALISM 



A. Supersymmetric quantum mechanics 



In this paper we consider SUSY in the one-dimensional 



KGE (1.4), and especially in the corresponding eigen- 
value problem 



where 



Hcj)j(x) — ajj(/)j(x) 



H = -dl + V(x) 



(2.1) 



(2.2) 



The boundary conditions will be specified later. In so 
far as the interest centers on the time-independent prob- 
lem (2.1) and the spectrum, the Schrodinger equation, 
to which reference is usually made , is included if we 
simply relabel uj 2 i— » E. 

If there exists another system 



H = -dl 



V(x) 



(2.3) 



with the same spectrum (or the same spectrum apart 
from one state), and moreover if the states in the two 
systems are related by 



where 



4>{x) = Acj)(x) 

A = d x + W{x) 
-A 1 =d x - W(x) 



(2.4) 



(2.5) 



then the two systems are said to be SUSY partners. In 
particular, if <fij(x) is an eigenfunction of H, then <fij(x) 
(provided it does not vanish) is an eigenfunction of H 
with the same eigenvalue. Normalization is deferred to 



Section III 



In order for (2.4) to preserve the spectrum, one needs 



AH = HA , 

from which it follows that 

V{x) = W{xf - W'{x) + Q 2 
V{x) = W{xf + W'{x) + Q 2 



(2.6) 



(2.7) 



with W{x) (called the SUSY potential) as in and for 
some constant ft 2 . Since both V and V have to be real, 
W and f2 2 are also real. Moreover, the Hamiltonians can 
be represented as 



H = A^A + n 2 
H = AA^ + fl 2 



(2.8) 



The two partner systems can be put into one linear 
space by introducing Pauli spinors, with H and H asso- 
ciated with 1 ± er 2 and A, A^ associated with a ^ 



Upon reversing the sign of W, (a) V «-> V [cf. (2.7)], 
and (b) A <-> — A^ [cf. ( |2.5| )]; thus the mapping from H 
back to H is (up to a sign) achieved by A* . Note however 
that the mapping is the "inverse" only in a loose sense: 
A^A is not the identity but H - Q 2 [cf. (E 



We may regard (2.7) as a Riccati equation for W in 
terms of the given V. For |x| > a, both V and V vanish, 
so W 2 = — O 2 . The first-order Riccati equation for W 
can satisfy two boundary conditions (at x = ±a) only at 
special values of f2 2 ; this condition becomes familiar if 
we define a generator <&(a;) by 



W(x) 



(2.9) 



(2.10) 



Then (f2/7|) implies 

H<f>(x) = n 2 $(x) . 

B. Boundary conditions 



All the above may be regarded as a review of the fa- 
miliar SUSY formalism for the Schrodinger equation [^JJ] 
if we put E — uj 2 , Eq = fl 2 and in particular shift V(x) 
so that Eq — 0. Conventionally the discussion refers 
to wavefunctions which vanish at infinity (or, more pre- 
cisely, are square- integrable). Here we consider all eigen- 
functions in T, including both NMs and QNMs, with the 
former in the upper and the latter in the lower half of 
the frequency plane. 

We should check immediately that <fi £ T implies G T. 
For x > a, if <j>{x) = Ce luJX then 0(a) = A(f>(x) = 
(iuj + W+)Ce MX , where 



W{x=±a) = W± 



(2.11) 



are the constant values for x > a and x < — a respec- 
tively. Thus and always satisfy the same type of 
boundary conditions, and the number of eigen stat es in 
T is preserved under SUSY: n(oj) — h(oj) [cf. ( |l.2| )1 - 
except when A or /U destroys a state, to be discussed 
below. 



C. Generator 



The various SUSY transformations are related, in a 



one-to-one manner, to solutions of (2.1C) for the gener 



ator $. First, suppose fl > 0, so that O is real. Then 
outside /, $ is oscillatory: either complex (e.g., e inx ), 
inadmissible since it leads to a complex W; or real (e.g., 
sin fix), inadmissible since its nodes lead to singularities 
in W. Thus, fl 2 < 0, and we denote K = 

At each spatial extreme (\x\ > a), $ is in general a 
sum of increasing and decreasing functions, i.e., 
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ce 



K\x\ 



de 



-K\x\ 



(2.12) 



If both c, d 7^ (to be called the mixed type), then the 
logarithmic derivative is (e.g., for x > a) 



W{x) = —K - 



2dK 



-2Kx 



(2.13) 



so that V — V + 2W' acquires an exponential tail. (In 
the special case O = 0, the tail is not exponential but 
asymptotically inverse-square.) Thus, if we insist that 

V also has finite support, the mixed type is not allowed 
and at each extreme $ must be either purely decreas- 
ing (<I> oc e~ K)[X \ denoted as D) or purely increasing 
($ oc e K)[X \ denoted as I). Outside /, the logarithmic 
derivative is then exactly ±K, so W = 0, implying 

V = 0. When both extremes are considered together, 
$ must be one of three types, conveniently labelled with 
the parameter 



X 



sign(W+) — sign(VF_) 



(2.14) 



where x = +1) — 1; an d respectively for the DD, II, 
and DI/ID cases (in obvious notation). In the DI/ID 
case, the generator is purely incoming from one ex- 
treme and purely outgoing to the other, hence is a total- 
transmission mode (TTM). The relaxation to allow ex- 
ponential tails for V and/or V will be briefly mentioned 



in Section VI D 



The gen erato r is annihilated by SUSY: A$ — 0, triv- 
ially from ( |2.5| ) and (|1]). Furthermore, since reversing 
the sign of W interchanges the partners [cf. below (2.8)], 
in view of (2.9) the transformation from H back to H 
is generated by $ — this is guaranteed to be an 

eigenfunction of H, also with eigenvalue il 2 . (Despite 
the notation, $ is not the SUSY partner of <l ^ A3>.) 
The boundary conditions for $ = $ _1 interchange D 
and I, so the reverse transformation is characterized by 

X = -x- 

These observations allow a simple statement of the 
changes in the number of states when H i— > H. If x — 1 
[cf. below ( pli| )], an NM $ is destroyed at Q = iK 
[A{iK) = 1] and a QNM <l is created at -Q = -iK 
[A(-iK) = -1]. If x = -1, a QNM $ is destroyed at 
fl = -iK [A(-iK) = 1] and an NM $ is created at 
-Q = iK \A{iK) = -1]. If x = 0, no eigenstates of V 
are created or destroyed [A(iK) — A{—iK) — 0], since 
$ and $ are TTMs rather than NMs or QNMs. Thus, 
all three cases satisfy 



A(iK) = -A(-iK) = x 



(2.15) 



where we emphasize the convention 57 = ±iK with 
K > 0. The cases x = ±1 lead to Hamiltonians whose 
spectra in T differ by one state (said to be essentially 
isospectral), whereas the case x = leads to Hamiltoni- 
ans whose spectra in T are identical (said to be strictly 
isospectral). 



These remarks provide a more complete picture of the 
mapping of eigenstates in T under SUSY: states do not 
simply appear/disappear, but are mapped to the mirror 
point in the complex plane. 

Not all NMs, QNMs or TTMs are eligible as the gen- 
erator. First, fl 2 must be real, which restricts QNMs to 
zero modes. Secondly, $ cannot have nodes, or else W 
would acquire singularities. In the case of NMs, this re- 
stricts $ to the ground state. For QNMs, nodes are not 
required by general theorems. At least for repulsive po- 
tentials, each eigenfunction can have at most one node or 
antinode; thus, for a symmetric repulsive V, even-parity 
eigenfunctions can have no nodes. There is consequently 
much more freedom in choosing a QNM (as opposed to 
an NM) as the generator. Some general statements con- 
cerning nodes in QNMs are given in Appendix |X| 

III. ORTHONORMALITY 

A. Orthonormality for NMs 

For conservative systems, SUSY preserves orthonor- 
mality. There are two issues: orthogonality is preserved 
because the transformed NMs are eigenvectors of the self- 
adjoint operator H ; and normalization is preserved if the 
transformation is changed to 



with 



<f>j >-> h = N ih = N i A< t>3 



N -2 = fefe) 

3 mh) 



(3.1) 



(3.2) 



Eq. ( p.l| ) applies to each eigenstate j other than the gen- 
erator $ itself, namely the ground state. It is readily 
shown that N~' 2 — u 2 — O 2 , a result that can also be read 
off as a special case of the derivation below for states in 



T. Since w 2 -S! 2 is the eigenvalue of A^A [see (|j|], ( |3T| ) 
can be written in the operator form 



l = A(A^A) 1/2 <p 



(3.3) 



valid for any state <f>, not just frequency eigenfunctions. 
(This and similar formulas below are restricted to the 
subspace orthogonal to $.) This makes it formally easy 
to verify the preservation of inner products: 



(^) = (^I(AU)- 1 / 2 At A (^ A)-^ 2 \ 



(3.4) 



However, when operating on a general wavefunction <fi, 
the factor (At A) -1 / 2 can only be evaluated by project- 
ing <f> onto the eigenfunctions, and scaling each compo- 
nent by Nj. Thus, in practice, the significant result is 
the evaluation of this factor. We now generalize these 
concepts to states in T, in particular QNMs. 
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B. Normalization and inner product for QNMs 

It is necessary to digress and review the concepts of 
orthogonality and normalization for QNMs. The central 
issue is that with the OWC, H is not self-adjoint in the 
usual sense, and different QNMs are not orthogonal un- 
der the usual inner product. Likewise, the norm J \<f>j \ 2 dx 
is divergent, since the wavefunction grows exponentially 
at infinity. 

An appropriate normalizing factor for QNMs was first 
introduced by Zeldovich and later generalized and 
applied to other situations including models of lin- 
earized waves propagating on a Schwarzschild back- 
ground p5[: 



(3.5) 



bj ,cj)j) = 2tUj / <fij (x) 2 dx 

J —a 

+ *[<M-a) 2 + 0» 2 ] 



This expression goes as cf> 2 rather than \4>j\ 2 , and is in 
general not real. The limits of the integral and the sur- 
face terms can be shifted from ±et to any b±, where 
±b± > a, without affecting the value of ( |3.5| ). This defi- 
nition also applies to NMs: the surface terms vanish if we 
take b± — > ±00, recovering the conventional norm apart 



from a factor of 2luj. In the QNM case, (3.5) is the cor- 



rect normalizing factor in the sense that, e.g., under a 
perturbation V ^ V + AV, the complex QNM eigenval- 
ues change by 



AK 2 ) 



Jcj) j (x) 2 AV(x)dx 



(3.6) 



Since one no longer has positivity, there is the possibility 
that (<pj ,4>j) = 0. This exceptional case can be separately 
taken care of 26 1 , and some interesting aspects are dealt 
with in Section IV. 



To go beyond the normalizing factor and discuss an 
analog of orthogonality, one has to first regard each state 
as a two-component vector i/> = {ifr 1 ,ifP , ) T = c?t"0) T , 
which is most easily motivated by noticing that the dy- 
namics requires two sets of initial data. The space T is 
then defined as all ip satisfying 



^ 2 (±a) =Td x ^(±a) 



(3.7) 



For an eigenstate, <pj — (<f>j,—iujj(f)j) T , and (1.7) would 
follow from (3.7). 

In terms of the two-component vector, one can define 
a bilinear map @,M 



(V, <f>) = i y [V- 1 (») + ^ 2 (x)^ 1 (x)] dx 

+ ^ 1 (-a)^ 1 (-a)+V^ 1 (a)0 1 (a)] } , (3.8) 



to take the place of the usual inner product. For an eigen- 
function, (<pj^<pj) agrees with (3.5). The dynamics can 
be written in the first-order form idt4> = with 



H 







1 



dl-V 



Importantly, Ti. is symmetric: 
(W^,0) = (^,W^) 



(3.9) 



(3.10) 



in the proof of which the surface terms generated upon 
integration by parts exactly cancel against those in ( |3.8] ). 
The relation ( 3.10| ) is the analog of self-adjointness, and 
leads to the usual proof that for two eigenvectors, 



{4>k,4>j) 



(3.11) 



whenever ujk 7^ Wj. Provided that (<f>j,cj)j) ^ Pq| . one 
can normalize thes e eigenfunctions in the usual way, i.e., 
by requiring ( 3.11 ) to be 2wk&kj in general [cf. (3__5) for 
this factor]. We hence forth refer to this property as or- 
thonormality [and to ( 3.11 ) alone as orthogonality]. It 
also follows trivially that, provided this orthonormal sys- 
tem is complete (which is the case under fairly broad 
assumptions; see Section |lB[ ), time evolution is given by 



cj)(t) = J2 a j(f>je' 



(3.12) 



generalizing (1.8) to two components, and 

The preservation of orthonormality under SUSY should 
therefore be sought in terms of the bilinear map (3 



(3.13) 



C. Normalized SUSY transformation for QNMs 

We first present a derivation of orthonormality that 
does not explicitly require the two-component formal- 
ism. With orthogonality already guaranteed by (3.11), 
it remains to compute the normalizing factor 



(3.14) 



6 j ,</> i ) = 2£J i / [{d x + W)<PjY dx 

J —a 

+ i 4>j{-a) 2 + <^(a) 2 l 



Integrate by parts to convert (d x (j)) 2 to —(d 2 (f))(f) plus 
a surface term, express the second derivative in terms 
of V — u> 2 by means of the eigenvalue equation, and 
write the potential as V = W 2 - W + n 2 . Then, 
apart from a term oc uj 2 — f2 2 , the integrand becomes 
a total derivative d x {W4> 2 ). Using lU(±a) 2 = — f2 2 and 
d x <pj{±a) — ±iu!j4>j(±a) then leads to 
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(3.15) 



Incidentally, the conservative case (nodal conditions at 
the ends of the interval [M ) is recovered by simply drop- 
ping all surface terms. 

Since the ratio (3.15) is the eigenvalue of A' A, we can 
again write the normalized transformation for each eigen- 
f unction as (3.3). 



D. SUSY for two-component form 

For eigenstates, the two components are trivially re- 
lated, but in order to perform SUSY transformations on 
general wavefunctions in T (e.g., given a time-dependent 
state, to find its partner at all times), the second compo- 
nent must be considered explicitly. 

Since SUSY must commute with time evolution and 
4> 2 = dt4> , both components must transform in the same 
way. Thus, the (unnormalized) SUSY transformation on 
two-component vectors is A = diag (A, A), which satis- 
fies 



(„4V» = W^V) 



(3.16) 



where A 1 = diag {A\ A^). j n deriving (3.16), one has to 
integrate by parts: d x changes sign so that A turns into 
AJ . The surface terms are seen to work out by using, 
e.g., 4> 2 (a) = dt<p 1 {a) = — d x cj) 1 (a), and the known values 
ofW±. Note that A* A = (H-Q?)l, AA^ = {H~n 2 )t, 
i.e., the products do not relate to the two-component Ti.. 

With (§T6|), it is straightforward to show that the nor- 
malized SUSY transformation 



0^ j> = A(A^A)- 1/2 cf> 



(3.17) 



defined on the subspace orthogonal to <I>, prese rves the 
bilinear map, in a manner that exactly parallels (3^4). In 
the exceptional case of SUSYs that generate a doubled 
state (see Section |i"v| ), the subspace has to exclude the 
two states on which H — il 2 vanishes. 

The linear-space structure for open systems (e.g., the 
replacement of inner products by bilinear maps) has an 
intrinsic geometric meaning for all outgoing states, not 
just QNMs [^7j . It is therefore pleasing that this struc- 
ture is preserved by SUSY, a superficially unrelated con- 
cept. 



IV. JORDAN BLOCKS 

A key concept in SUSY is the preservation of the spec- 
trum. However, dissipative systems (such as waves satis- 
fying the OWC) admit a spectral property not found for 
conservative systems. In terms of the Wronskian J(uj) 
to be defined below, this is exhibited as an Mth-order 



zero (M > 1). Such a multiple zero emerges naturally 
when M QNM eigenvalues coalesce as system parame- 
ters are tuned, so that M—l eigenvectors are "lost" pi| 
and must be replaced by other degrees of freedom. Thus, 
the Hamiltonian cannot be written as a diagonal matrix 
in the (biorthogonal) basis of eigenstates, but can only 
be decomposed into (Jordan) blocks |l(| of size M x M. 
When this happens, ((j)j,<pj) will vanish for some j, in- 
validating the formalism in Section III [see, e.g., (|3 



J 

and (3.13)]. These issues have been discussed in detail 
with reference to waves in open systems p3]. 

The simplest example of a JB (with M = 2) is a har- 
monic oscillator going through critical damping. The 
eigenvalues ui± — ±wr — coalesce when cjr — > 0. With 
one eigenvalue lost, the dynamics is not given by a sum 
of eigenfunctions with time dependence exp(— iui±t), but 
by only one such eigenfunction, plus another term whose 
time dependence carries a prefactor t. 

Our purpose in this Section is to establish that SUSY 
maps a JB in H into a JB in H, preserving the order M 
except when the eigenvalue coincides with ±fi. 



A. Wronskian 

In this subsection we introduce the Wronskian J(u>), 
define JBs in terms of its multiple zeros and describe the 
mapping of JBs under SUSY by a relation between J(uj) 
and its counterpart J{ui). 

In the original system H , define solutions of the wave 
equation f(u>, x) and g(u>, x) satisfying the boundary con- 
ditions /(oj, — a) — 1, f'{u),—a) — —ilu, g(ui,a) = 
1, g'(oj,a) = iu>, where ' = d x . The function values 
are arbitrary normalizations, while the derivatives im- 
pose the OWC on the left and right respectively. An 
eigenstate <j>j in T satisfies the boundary condition on 
both the left (as for /) and the right (as for g): <fij cx 
f(u>j,x) oc g(u>j,x). Thus, the zeros of the (position- 
independent) Wronskian 



J M = /V, x)g{u>, x) - f(u, x)g'{u>, x) 



(4.1) 



identify the eigenvalues in T. It can be shown [see (4.7) 
below] that (jfij, <pj) oc dJ(u>j)/du>, so an Mth-order zero 
of J [M > 1) corresponds to the generalized norm being 
zero, and is precisely the JB phenomenon that we wish 
to investigate. 

It is natural to generalize the definition of n(uj) to be 
the order of the zero, viz. 



1 £ dJ (lo' ) / du)' , 



(4.2) 



on a contour of winding number +1 enclosing u>. This 
definition makes it clear that the total number of states 
(but not necessarily of eigenstates) within a contour is 
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preserved under continuous changes of the system pa- 
rameters. We note for future reference (see Section VI D| ) 
that poles of J (which can only occur if V does not have 
finite support) count as negative values of n. 

Now consider the analogous construction in the part- 
ner system H, obtained for example by using an NM $ of 
H as the generator, i.e., for \ = 1. By our convention, $ 
is associated with a frequency SI = iK, and W± — =FzSl. 
The SUSY transformation gives pi 



f(w,x) = (d x + W)f(u,x) 
g(w, x) = (d x + W)g(u, x) , 

leading to the Wronskian 

= /'(w, x)g(uj, x) - f(u, x)g'{uj, x) 



(4.3) 



(4.4) 



When (4.4) is written out using ([4.3]), some terms can- 
cel by using J' = 0, the second derivatives can be elim- 
inated by the defining equation, and V is expressed in 
terms of W and SI 2 . Some arithmetic then leads to 



j u (w) = {u? - si 2 ) j(u) 



(4.5) 



This Wronskian is however unnormalized (as indicated by 
the superscript), since only Cf and Dg satisfy the nor- 
malization conventions at —a and +a respectively, where 
C = —D = i(oj — SI) -1 . Thus the normalized Wronskian 
J(w) = CDJ u (u) is 



J(w) 



n 



(4.6) 



The central result (4.6) neatly summarizes the corre- 
spondence between the two spectra for x = 1. [Similar 
formulas for the other cases can all be consolidated by 
changing SI t— > ixK in (4.6), and will not be separately 
discussed.] (a) For uj ^ ±fi, the spectra of H and H 
are the same: a JB of order M in H maps to a JB also 
of order M in H at the same frequency. This should 
be no surprise since up to the point of coalescence (for 
M > 2), the eigenvalues of the two systems are guaran- 
teed (cf. Section [n]) to be in one-to-one correspondence; 
in a sense, the result here merely demonstrates that the 
limit is not singular, (b) Moreover, at the special fre- 
quencies ±S1, M states at SI (of which only one is an 
eigenstate) are mapped into M—l states at SI plus one 
state at — SI. In other words, we recover (2.15) even for 
JBs, i.e., even when n, n > 1. Anticipating the possibil- 
ity of poles in J (cf. Section VI D), we note that (4.6) 
implies that such poles are also preserved under SUSY, 



and would be accommodated by (1.2), (1.3), and (2.15) 
with negative n, n. 

Incidentally, (3.15) on the change in normalization un- 
der SUSY follows simply from (4.5), since the bilinear 
map is related to the Wronskian by p4],pfj| 



dJ(u) 



did 



(4.7) 



B. Doubling of states by SUSY 



We have seen that (say for \ = 1) one has n(— iK) — 
n(—iK) = 1. Where this increases from to 1, the sit- 
uation is straightforward — a QNM is created. When 
the increase is from 1 to 2, the situation is more sub- 
tle and merits a detailed examination. [The general case 
where this increases from M to M+l (M > 1) will not 
be shown.] 

Consider a system H with an NM $ at SI = iK and ac- 
cidentally also a QNM ^ at -fl' « -£l. If SI' ^ fi, there 
are two corresponding QNMs in the -ff-system, namely 
$ = at -SI and ^ = A$j at -fi'. Now tune the 
parameters of H so that SI' — > SI; in the limit we must 
have i&j cx $ [ p6p§| ], as is readily verified. The propor- 
tionality constant can be evaluated by 



*i = ttj(-a) 
$ $(-a) 

' *(<*) 



= 2ifi# J -(-a)$(-a) 
-2ifi*j(a)$(a) . 



(4.8) 



The agreement of these two expressions can also be seen 
without invoking SUSY. One notes that €> and SPj, be- 
ing eigenfunctions of TL with distinct eigenvalues, are or- 
thogonal. In the bilinear map, the integral term van- 
ishes because the frequencies are opposite, leaving only 
the surface terms. Thus one finds = (*&j,$>) = 
i[^-(-a)$(-a) + %(a)$(a)]. 

Now the frequency — SI in H must be associated with 
a doubled state. This can be seen in two ways: (a) until 
the limit SI' — SI, the re are two distinct states; (b) from 
the key relation (4.6), J has a second-order zero. With 
the two QNMs collapsed into one, there has to be another 
basis vector, to which we now turn. 

Using the normalization of /, one has ^j(x) — 
®j(-a)f{-n,x), implying Vfe) = ^(-^/(-fi, x). 
[Analogous formulas for <E> and $ follow from (|4.8|) .] The 
double zero of J(to) at u> = —SI means that f(u),x) sat- 
isfies the OWC at x — a not only at u> — —SI, but also 
to first order away from the zero. This makes it plau- 
sible that, in the QNM expansion, d w f(u>, x)\_q takes 
the place of the "missing" eigenfunction when $ and ^fj 
coincide, which has been confirmed in detail |^6|. One 
thus defines, for an arbitrary a, a pair of functions *&j, n , 
where n = 0, 1 is an intra-block index: 

*J,o(aO = 

¥ j , 1 {x)=* j (-a)d u f(w,x)\-a+a9 j (x) , (4.9) 
and the second function satisfies 
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(4.10) 



V. EXAMPLES 



Using this, one verifies that an outgoing solution is given 
by *j,i(t) = (%i - it^j t o)e int . This time dependence 
shows that the associated second component should be 
'J 2 '• x = — i{— 0*j,i + ^ jfl)- The prefactor t and its coun- 
terpart if any in the ii-system will be further explored 
below. 

Next consider the bilinear map and normalization for 
these functions. For a double zero one always has 
(*j,o, *j,o) = 0; cf. (O) [a choice of a in (O) such that 



also = would be useful in wavefunction ex- 

pansions]. The JB is normalized by one overall factor, 
the bilinear map between the two basis states |2q] : 



~^d 2 J u (uj)/duj 2 
—dJ{uj)/duj 



J -n 



2Q 



(4.11) 



where in the numerator we have used the result analogous 



to (4.7) for a double zero. 

Finally the reverse transform generated by — |3l| l 
satisfies the following properties, (a) A'^jfi cx = 0. 

(b) Hence in A^jift), the term cx te int is annihilated. 

(c) The remaining term in A^j t i is d&j, readily seen by 
observing that 

(H - n 2 )(A^ jt i) =A\H- 2 )* 3 ,i = A f (-2iW jfi ) 
= , (4.12) 

so that A^j i is an eigcnfunction of H with eigen- 
value O 2 . In particular, the time dependence is e 
without any prefactor t. A straightforward computation 
shows that c = — 20. 

Since A^A is not the identity, two different states in 
the ii-system can be associated with \F,,i. The first is 
A^j t i as discussed above. The second is the SUSY pre- 
image of under A, which is readily found by noticing 
that 

% i(x) = tf^-a) x d u [(d x + W)f(u,x)]- a 

= J A* J (-o)5L,/(w,x)|_ n . (4.13) 

However, ^fj(— a)<9 w /(u;, x)|_n is not outgoing, since 
J{oj) only has a first-order zero at — fi, and consequently 
f(uj,x) satisfies the OWC at x = +a only at —CI, but 
not to first order away from it. 

These remarks completely resolve the puzzle related 
to the prefactor t in the time evolution in the ii-system. 
Namely, of the two corresponding wavefunctions of H, 
one has an exponential time dependence, while the other 
is not in T. 



The formalism developed in this paper is general, in 
that given V(x) with finite support, SUSY partners can 
be constructed whenever there are nodeless generator 
candidates. Nevertheless, some simple examples will suf- 
fice for an illustration. 

Let V be a square barrier of height Vq on i; without 
loss of generality henceforth a = 1. In the even sector 
there are two zero modes for small values of Vo; e.g., for 
Vq = 0.16 they occur at cji = — 0.18K, uj 2 — — 2.500i. 
The wavefunctions are <j)j(x) = coshojir within / (with 
ai = 0.242, «2 = 2.506), and a real exponential for 
| or | > a; clearly each <pj has no nodes. 

This example already illustrates, as unexceptional, the 
existence of several nodeless zero modes — any one of 
which can be used as the generator <E>. We choose the 
state at O = u>\ ; since $ is purely increasing at both ex- 
tremes, x = — 1- Figure 1(a) shows V (solid line) and V 
(dotted line). Figure 1(b) shows the spectra in the com- 
plex w-plane: most eigenvalues are common (crosses); 
one QNM exists only in H (circle), while one NM exists 
only in H (triangle) — the two systems are essentially 
isospectral. Another essentially isospectral partner can 
be constructed using the state at u)2 as the generator. 
The reverse transformations are characterized by x = 1. 

As Vo increases, the two zero modes move closer to- 
gether and merge at Vo = 0.291, forming a JB in H (with 
M = 2). This example illustrates the general results of 
Section |y| (with H and H interchanged). 

As another example, consider a symmetric multi-step 
square barrier, with V(x) being —10 for |x| < 0.1, 1 
for 0.1 < |x| < 1 and zero for |x| > 1 [solid line 
in Figure 2(a)]. In this example, there is a TTM at 
Q = — 0.990i [square in Figure 2(b)] — in fact a par- 
ity doublet with one propagating from the left (TTMl) 
and one propagating from the right (TTMr ). We choose 
the former as the generator for a x = transformation. 
The partner potential V is shown by the broken line in 
Figure 2(a). Since the generator is not symmetric, nei- 
ther is V. The states in T are exactly preserved [crosses 
in Figure 2(b)]. 

In this example it is interesting to consider not just 
the states in T, but also TTMs (see also Section VIC). 
By arguments similar to those in Section O, one TTMl 
$ is destroyed and one TTMr $ = is created at 
the same frequency O [square in Figure 2(b)]. However, 
in this example, because V is symmetric, there is also a 
TTM R l'(x) = $(-x) at O, and its partner <l = A^ is 
again a TTMr in the ii-system. Thus, in the ii-system, 
there is a doubled TTMr state at O. The situation can 
again be analysed in terms of the double zero of a Wron- 
skian Jr(w), but in this case Jr and Jr refer to TTMr's. 
For example, Jr is now defined in terms of a function / 
satisfying the outgoing-wave condition at x = —a and 
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a function g satisfying the incoming-wave condition at 
x = a. The obvious adaptation of the discussion in Sec- 
tion |[V] will not be given. 



VI. DISCUSSION 

A. Summary 

In this paper, we have extended the usual discussion 
of SUSY as a relation between NMs of partner systems 
to include the QNMs as well. By viewing all these to- 
gether in the space T, a more complete picture emerges. 
For example, in the usual discussion for NMs only, essen- 
tially isospectral transformations are said to lose or gain 
one state; now we see that (for \ — ±1) when an NM 
appears (disappears), a corresponding QNM disappears 
(appears). 

Furthermore, we have shown that the nontrivial linear- 
space structure for QNMs and any possible JBs are pre- 
served by SUSY — the latter being a feature not found 
in conservative systems. 

QNMs differ from NMs in two further regards. First, 
they have complex frequencies; nevertheless, even with 
twice as many constraints, matching the spectra turns 
out to be not any more difficult. Second, QNMs need 
not have an increasing number of nodes, and it is often 
possible to find several nodeless QNMs which generate 
distinct SUSY transformations — whereas the analogous 
operation for NMs would restrict the generator to the 
nodeless ground state. 

These wider perspectives are gained only because at- 
tention is paid to the Klein-Gordon rather than the 
Schrodinger equation, since the concept of outgoing 
waves has no meaning in an equation that is first-order 
in time. 

Two further important properties are also preserved, 
(a) If V has a singularity say at x = ±a (e.g., a step), 
then V will have the same type of singularity, but with 



opposite sign, as can be seen from (2.7) by noticing that 
the most singular part is W. (b) If V has finite support, 
then provided $ is not of the mixed type, V would like- 
wise have no tail. These two properties are precisely the 
conditions for the eigenstates in T to be complete 
Thus, SUSY maps a complete basis to a complete basis 
(if, for the \ = ±1 cases, $ t—> $ is included as well). 



B. Inversion 

This work partially answers the question of QNM in- 
version. It is well known ]33[ ] that on a finite interval, 
two sets of real NM frequencies uniquely determine the 
potential V. Does one set of complex eigenfrequencies 
in r also uniquely determine VI The answer is negative: 
there can be strictly isospectral potentials if a TTM with 



purely imaginary frequency exists and can be used as a 
X — generator; Figure 2 provides one such example. 
However, the following scenario is not yet ruled out. If 
we consider a half-line problem < x < oo (say corre- 
sponding to the radial variable in a 3-d system) , imposing 
a nodal condition at x = and the OWC for x > a, can 
one set of QNM frequencies uniquely determine the po- 



tential? SUSY transformations (2.4) do not directly re- 
solve this possibility — for which there is some numerical 
evidence |34j — since these one-sided systems do not fea- 
ture an analog of TTMs with which one could construct 



strictly isospectral partners. Moreover, by (2.7) and fl2.9| ) 
the nodal condition maps a regular V to V ~ 2/x 2 for 
x — > + (generalizing the well-known result that SUSY 
increases the angular momentum by one unit in the hy- 
drogen atom). It would therefore be interesting to see 
if an enlarged class of transformations can address this 
question. 



C. Total-transmission modes 

The present paper refers in the main to states in T, i.e., 
states that satisfy the OWC at both extremes. One could 
also develop the same formalism for TTMs; see, e.g., the 
end of Section fv]. Note that we here consider TTMs as 
states on which SUSY acts, rather than as (x = 0) gen- 
erators. To be more specific, the SUSY transformation 
acts on the space Tl of TTMl's 



<j>(x) 



oo 



(6.1) 



or (equivalently, under lu t— > —uj) the space Tr of TTMr's 



<f>(x) 



(6.2) 



One significant difference is that a \ = ±1 transforma- 
tion preserves the eigenvalues in these spaces, whereas a 
X = transformation may shift one state. 



D. Tails 

So far we have only considered potentials with finite 
support. We end with an outline of some issues that 
arise for long-range potentials. 

The most obvious difference is that if V and/or V are 
allowed to have tails that decay exponentially or slower 
in \x\, then the mixed- type SUSY transformation is al- 
lowed, and there is a continuous choice of generators. 

However, there are a range of more subtle issues. First, 
the very definition of an outgoing wave requires care. At 
the numerical level, special treatment is needed to ensure 
convergence pa], but there are matters of principle as 
well when the OWC <f>(w, x) ~ e + lu \ x \ is imposed only as 
\x\ — > oo. In that limit, the condition as stated becomes 
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vacuous for Im uj < 0, since admixture of an (exponen- 
tially smaller) incoming solution would not alter this be- 
havior. Rather, it is necessary to define the OWC in the 
upper half-plane in uj and analytically continue to the 
lower half-plane. A wave 4>{uj,x) is incoming if uj,x) 
is outgoing this is equivalent to saying that incom- 
ing waves are defined first in the lower and continued to 
the upper half-plane. The necessity for these procedures 
makes it possible that (at certain singular points uj) a 
wavefunction can be both incoming and outgoing. 

Singularities in the one-sided functions ]36| ] can only 
occur because of the need to integrate the defining equa- 
tion over an infinite range. Those values of lu for which 
these functions (say the left function /) are singular are 
said to be anomalous; if the potential is not oscillating, 
anomalous points can only occur on the imaginary axis 
- precisely where possible SUSY generators are to be 
found. The case of an exponential tail V(x) ~ V\e~ x \ x \ 
is of particular interest because the anomalous points 
(at uj m = —imX/2, m = 1,2,...) can be studied by the 
Born approximation, which in this case is equivalent to a 

More generally, if 



power-series expansion in z 



V(x) = v ke 

k=l 



- A I x 



k\\x\ 



(6.3) 



then for particular choices of the coefficients Vk , one (or 
more) of the generically anomalous points uj m may turn 
out to be regular — a situation wc refer to as miraculous. 
The anomalous and miraculous properties for the one- 
sided functions / and g are inherited by their Wronskian 
J, which is central to the formalism. These concepts have 
been discussed in relation to a particular application |37| . 

As far as SUSY is concerned, we only make one re- 
mark: such singularities lead to poles in J(w), thus are 
associated with negati ve val ues of n(u>), and are related 
in SUSY partners by fl2 15|) . The Poschl-Teller poten- 
tial V(x) — V sech 2 (a:/b) |38| illustrates many interesting 
features, including the existence of double poles of J{uj) 
which can merge with two zeros and lead to "nothing" - 
a potential that has total transmission at all energies |39) . 

Anomalous points are exceptional (in the case of ex- 
ponential tails being a discrete set of measure zero in the 
w-plane) and miraculous cancellation of singularities dou- 
bly exceptional. One might therefore think that these are 
not important. Surprisingly however, the problem of lin- 
earized gravitational waves propagating on a black-hole 
background jh| is precisely miraculous in this sense, at 
the so-called algebraically special frequency fl El[. A 
generator at f2 leads to a SUSY transformation that ex- 
actly relates the axial and polar sectors, which are there- 
fore isospectral in T. Among the more intriguing results 
is the following p7\ : in the polar sector there is a mode 
at £1 that is simultaneously a QNM and a TTM [i.e., 
at radial infinity it is purely outgoing but into the event 
horizon it is both outgoing and incoming] , while no modes 



exist in the axial sector. 

The subtle and perhaps counter-intuitive nature of 
these concepts demands a separate and rigorous exami- 
nation, to which the foregoing is meant only as a preview 
- and as further illustration of the utility of SUSY in 
dealing with waves in open systems. 
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APPENDIX A: NODES IN QNM 
WAVEFUNCTIONS 

Nodeless eigenstates play a special role in SUSY: they 
are candidates for the generator $. It is therefore useful 
to highlight the differences between NMs and QNMs in 
this regard, especially to contrast with the well-known 
property that there can be only one nodeless NM. 

First of all, we show that for QNMs with Re uj ^ 0, 
there can be at most one node or antinode. This is not 
surprising: since the eigenvalue is complex, the wave- 
function has a changing phase, and it would be "un- 
likely" that the real and imaginary parts (or their deriva- 
tives) would vanish together. To prove this formally, take 
the Schrodinger point of view, so that the eigenvalue is 
E = lo 2 with a nonzero imaginary part. Now consider a 
time-dependent QNM and suppose that it has nodes or 
antinodes at two points xi,X2- At these two points, the 
current 



J = t [<fi*(d x <P) - (d x cj>*)cj>) 



(Al) 



vanishes. Then, flux conservation implies that the to- 
tal probability in the interval [xi, X2] is constant in time. 
Yet the wavefunction is either growing or decaying, since 
Im E 0, which is therefore a contradiction. 

From the perspective of SUSY it is unfortunate that 
the above proof excludes the crucial imaginary axis. 
However, on that axis the statement remains valid for 
repulsive potentials, or more generally for potentials 
which are so weakly attractive that V — uj 2 is positive 
definite p2[ . Namely, let be a solution with two 
(anti)nodcs. By taking the real or imaginary part, we 
may assume <f> to be real. Now between two nodes <p 
would have an extremum, i.e. <f>"<f> < which is incom- 
patible with the KGE. Similarly, an antinode can only be 
a global maximum or minimum, precluding the presence 
of any other nodes or antinodes. 

Thus, except for the imaginary axis in the case of at- 
tractive potentials, QNMs can have at most one node or 
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antinode. For symmetric potentials, in the even sector 
x = is already an antinode, so there can be no nodes 
anywhere. 

For zero modes, i.e., QNMs with Re ui — 0, nodes are 
more "likely" : the eigenvalue is real and the wavefunc- 
tion has a constant phase (say purely real) , so a node re- 
quires only one condition, rather than two. Nevertheless, 
in contrast to the conservative case, the proof that there 
can be only one nodeless eigenstate can be bypassed. 

The interlacing nodal structure of NM eigenfunctions 
follows from well-established Sturm-Liouville theory. For 
the present purpose, we do not need the full apparatus. 
Consider, for simplicity, a finite interval [—a, a] and sup- 
pose there are two distinct nodeless eigenfunctions 4>i,<j>2, 
both real. Then they can both be chosen to be positive, 
which violates the orthogonality condition for NMs 



4>i(x)4>2{x) dx = 



(A2) 



We can attempt to transplant the argument to QNMs. 
For zero modes, the wavefunctions can again be chosen 
to be real, and if they are nodeless, positive definite. 



However, the analog of (A2) for two eigenfunctions with 
eigenvalues ojj = —ijj is 



(A3) 



-(71+72)/ 4>i (x) <j>2 (x) dx 

J —a 

+ [<M-a)<M-a) + ^(0)^2(0)] = 



Note in particular the signs of the two terms. With 
7j > 0, this condition does not preclude both eigenfunc- 
tions from being positive definite. 

Thus, we can make three remarks, (a) For NMs, there 
can be only one nodeless state, (b) For QNMs with 
Re u> = 0, there could be more than one state with no 
node, (c) For QNMs with Re w ^ 0, or with Rew = 
but V — oj 2 positive definite, each eigenfunction can have 
at most one node or antinode, and for symmetric poten- 
tials, every even eigenfunction is nodeless. 

Case (b) in particular opens up the possibility of mul- 
tiple SUSY transformations. 
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Fig. 1(a): A square-barrier potential V (solid line) and its SUSY part- 
ner potential V (broken line). Both potentials arc symmetric and only 
the x > part is shown. The SUSY transformation is constructed 
by using the state at SI - wi — — 0.18H [circle in Fig. 1(b)] as the 
generator. 
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Fig. 1(b): The complex oj-plane showing the QNMs common to both 
potentials (crosses); the mode present only in V (circle), which corre- 
sponds to the generator and the mode present only in V (triangle), 
which corresponds to <!> — <1>~ 1 . 
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Fig. 2(a): A multi-step potential V (solid line) and its SUSY partner 
potential V (broken line). The SUSY transformation is constructed by 
using the TTMl at O — — 0.990i [square in Fig. 2(b)] as the generator. 
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Fig. 2(b): The complex u;-planc showing the NM and QNMs com- 
mon to both potentials (crosses). The square indicates a TTMl and a 
TTM R in V, and a doubled TTM R in V. 
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FIGURE CAPTIONS 



Fig. 1 

(a) A square-barrier potential V (solid line) and its SUSY partner potential V (broken line). Both potentials 
are symmetric and only the x > part is shown. The SUSY transformation is constructed by using the state 
at O = Wi = — 0.18H [circle in Fig. 1(b)] as the generator. 

(b) The complex w-plane showing the QNMs common to both potentials (crosses); the mode present only in V 
(circle), which corresponds to the generator <!>; and the mode present only in V (triangle), which corresponds 

to <i = 

Fig. 2 

(a) A multi-step potential V (solid line) and its SUSY partner potential V (broken line). The SUSY transfor- 
mation is constructed by using the TTMl at ft = — 0.990i [square in Fig. 2(b)] as the generator. 

(b) The complex w-plane showing the NM and QNMs common to both potentials (crosses) . The square indicates 
a TTMl and a TTM R in V, and a doubled TTM R in V. 
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